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Abstract — In this paper some new experimental results about 
the statistical characterization of the non-line-of-sight (NLOS) 
bias affecting time-of-arrival (TOA) estimation in ultrawideband 
(UWB) wireless localization systems are illustrated. Then, these 
results are exploited to assess the performance of various 
maximum-likelihood (ML) based algorithms for joint TOA lo- 
calization and NLOS bias mitigation. Our numerical results 
evidence that the accuracy of all the considered algorithms 
is appreciably influenced by the LOS/NLOS conditions of the 
propagation environment. 

Index Terms — Ultra Wide Band (UWB), Radiolocalization, 
Bias mitigation. 

I. Introduction 

Wireless localization in harsh communication environments 
(e.g., in a building where wireless nodes are separated by 
concrete walls and other obstacles) can be appreciably affected 
by direct path attenuation and NLOS conditions. Various 
solutions for NLOS error mitigation in UWB environments are 
available in the technical literature 0, 0, 0, 0. A simple 
deterministic model, dubbed wall extra delay, is proposed 
in to estimate the bias introduced by walls. A non- 
parametric support vector machine is employed in |2| for joint 
bias mitigation and channel status detection; this approach 
exploits multiple features extracted from received signals in 
a non-statistical fashion. A few classification algorithms for 
LOS/NLOS detection are compared in [3|, where it is shown 
that the best solution is offered by a statistical strategy based 
on the joint probability density function (pdf) of the delay 
spread and the kurtosis extracted from the received signals. 
Finally, in [4] statistical models for the time of arrival (TOA), 
the received signal strength (RSS) and the root mean square 
delay spread (RDS) are developed and an iterative estimator 
for bias mitigation is devised. 

The contribution of this paper is twofold. In fact, first 
of all, the problem of joint statistical modeling of multiple 
features extracted from a database of waveforms acquired in 
a TOA-based localization system is investigated. Note that, 
as far as we know, in the technical literature only univariate 
models for bias mitigation have been proposed until now (e.g., 
see 0, 0). The use of multiple signal features in UWB 
localization systems has been investigated in [3] for channel 



state detection only and in 0, where, however, the considered 
features (namely, the kurtosis, the mean excess delay and the 
delay spread) are modelled as independent random variables. 
The second contribution is represented by a performance 
comparison of various maximum-likelihood (ML) estimators 
for TOA-based localization. In particular, unlike other papers 
(e.g., see 0, 0) we illustrate some numerical results referring 
to the accuracy of different localization strategies, rather than 
to the bias removal on a single radio link. 

The remaining part of this paper is organized as follows. 
In Section III] some information about our UWB experimental 
campaign and about the features extracted from the acquired 



data are provided. In Section III some estimation algorithms 
for UWB radiolocalization are described, whereas their perfor- 
mance is compared in Section [IV] Finally, some conclusions 
are given in Section [V] 

II. Experimental setup 
A. Measurement arrangement 

A measurement campaign has been conducted by our re- 
search group in 2010; it is worth mentioning that various 
databases providing a collection of sampled UWB waveforms 
acquired in experimental campaigns and useful for assessing 
the performance of localization algorithms are already avail- 
able (e.g., see 0, 0). However, our database has been specif- 
ically generated to assess the correlation between the NLOS 
bias error and various features extracted from the received 
signals, as it will become clearer in the next Paragraph. 

All the measured data were acquired by means of two FCC- 
compliant PulsON220 radios commercialised by TimeDomain 
and were collected in a database. Such devices are equipped 
with omnidirectional antennas, are characterized by a -10 dB 
bandwidth and a central frequency equal to 3.2 GHz and 
4.7 GHz, respectively, and perform two-way TOA ranging; 
they also allow to store the digitised received waveforms (a 
sampling frequency of 24.2 GHz and 14 bits per sample are 
used). Our measurement campaign consisted of two phases. 
First, the transmitter was placed in a given room (room A 
in Fig. [1]) and the receiver in an adjacent room (room B in 
Fig. [Til separated from the room A by a wall having thickness 
twaii — 32 cm (NLOS condition); in addition, the transmit 
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Figure 1 . Map of the environment of our experimental campaign. 



antenna was kept fixed, whereas the receive antenna was 
placed in A^° s = 174 distinct vertices of a dense square grid 
(the distance between a couple of nearest vertices was equal 
to 21 cm). In the second phase of our measurement campaign 
the transmitter and receiver were both placed in room B (LOS 
condition); then, the transmit antenna was kept fixed, whereas 
the receive antenna was moved on the same dense grid as 
in the first phase, to acquire the UWB signal in N%£L = 105 
distinct vertices. The distance between the two antennas varied 
between 1 m and about 5 m; larger distances were not taken 
into consideration since we were interested in indoor ranging 
only. It is important to note that the choice of the measurement 
scenarios described above is motivated by the fact that the 
UWB signals experienced similar propagation in both phases. 
For each grid position, besides the acquired waveform, a 
TOA estimate (evaluated by the PulsON220 devices) and the 
actual transmitter-receiver distance (evaluated by means of a 
metric tape with an accuracy better than 1 cm) were also stored 
in the database. The TOA estimates have been used to provide 
a common time frame to all the acquired waveforms; this has 
made possible the estimation of the mean excess delay and of 
other signal statistics in the signal processing phase. 

B. Statistical modeling of signal features 
The model 
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was adopted for the i-th link for both LOS and NLOS 
conditions. Here, n denotes the TOA estimated by the Pul- 
sON220 devices, based on energy thresholds and a go-back 
technique (see O for further details), cq is the speed of 
light, di denotes the distance between the transmitter and the 
receiver, b{ is the NLOS bias (in seconds) affecting the TOA 
measurement (the values taken on by this random parameter 
are always positive for NLOS links and null for LOS linksM) 
and Wi ~ A/"(0, o\ 4 ) is the measurement noise; in addition, 
the expression a 1 ^ { = 7<r T 2 jdf is adopted for the variance of 
the measurement noise, where 7 is a parameter depending 
on both the specific TOA estimator employed in the ranging 
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Table I 
Absolute value of the correlation coefficient between bi and 

EACH FEATURE OF THE SET {XiJ, j = 0, 1, ...., 5}. 



measurements and on various parameters of the physical layer, 
and j3 is the path-loss exponent (a known and fixed value is 
assumed for this parameter in both LOS and NLOS conditions 

my 

In the following we focus on the problem of estimating 
the bias 6, (affecting the TOA estimate Tj (fTh) from a set of 
Nf — 6 different "features" {xij, j — 0, 1, ..., 5} extracted 
from the received waveform r,(£). In particular, like in [2|, 
the following features have been evaluated for the set of the 
received waveforms: 

1) the maximum signal amplitude Xi t o — r ma x,i — 
max t \n(t)\; 

2) the mean excess delay Xi t \ — r m> , = /„ t 
parameter Si is defined below); 

3) the delay spread x ij2 = T dSti = / °°(t- T m ,i) 

4) the energy x if3 = e t = J °° \ri(t)\ 2 dt; 

5) the rise time Xi } 4 = t r ise,i = 
miri {t : |rj(f)|/max f \r. b (t)\ > 0.9} 
min{£: \n(t)\/ max t \n(t)\ > 0.1}; 

6) the kurtosis 2^5 = k% = -r T J T (\ri(t)\ — fi\ r \) dt, 

where p\ r \ = \ J T \n{t)\dt, of r] = ^ J T (\n(t)\ - 
p,\ r \) 2 dt and T denotes the observation time. 
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In the following the set {?"j(t)} of received waveforms is 
modelled as a random process, so that the above mentioned 
features form a set of correlated random variables; in addition, 
all of them are statistically correlated with the TOA bias. 
The last consideration is confirmed by the numerical results 
of Table |TJ which lists the absolute values of the correlation 
coefficients of the previously described features with the 
estimated TOA bias for both the LOS and the NLOS scenarios. 
From these results it can be easily inferred that not all the 
considered features are equally useful to estimate the TOA 
bias. For this reason and to simplify our statistical analysis, 
we restricted the set of features to {x i0 , Xn, 2^2} (Nf = 3), 
which collects the parameters exhibiting a strong correlation 
with the bias in the NLOS scenario. 

Note that, in principle, the bias bi is not influenced by the 
distance di, since it depends only on the thickness of the walls 
(or of other obstacles) encountered by the transmitted signal 
during its propagation; this is not true, however, for the above 
mentioned triple of signal features (see [1] for further details). 
Generally speaking, it is useful to derive a TOA bias estimator 
which is not influenced by the transmitter-receiver distance 
di. Therefore, in the attempt of removing the dependence 
of the features {x.^q, Xi t %, 2^2} on the link distance, we 



'in this case the probability density function (pdf) of 6j is f b (b) = 8(b). developed the models Xifl 
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Figure 2. Estimated joint pdf of the estimated TOA bias and the delay spread 
in NLOS conditions (above) and LOS conditions (below). 
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basis of our experimental results and, in particular, on the 
basis of the estimates of the joint pdf's {fb,x , fb,xi, fb,x 2 } 
referring to the three possible couples (pi, Xi j) with j = 0, 1, 
2 and i = 0, 1, ..., N acq - 1, where N acq = N^° s + N£ q ; here 
r max> T m an d T dl are random variables, whereas r^ s and r™ aa , 
are deterministic parameters having known values (all are 
assumed independent of di and thus also from the link index 
i). Given these models, the vector of distance-independent 
features x = [r ( } nax ,T™,T 7 Jl] can be evaluated from its 
distance-dependent counterparrl x = [r max , T m ,Td s } and can 
be used in place of it. However it turns out that parameters of 
x are less correlated with bi than those of x; for this reason, 
we decided to take into consideration both for localization 



purposes (see Sec. III-B2 1 



A complete statistical characterization of the estimated bias 
and of the 3 related signal features we consider is provided 
by the joint pdf fb.r max ,T m ,T dB (•) or, equivalently, by the joint 
pdf f br o T m T m (') (both can be estimated from the acquired 
data). Some marginalizations, like the one shown in Fig. [2] 
reveal interesting aspects: 1) a significant (limited) correlation 
between these parameters is found in the NLOS (LOS) case; 
2) the null region exhibited by the estimated pdf's is due 
to the fact that the TOA bias cannot take on values in the 
interval [0, t wa u/co], where t wa u is the thickness of the wall 
obstructing the direct path; 3) large values of the TOA bias are 
unlikely since they are associated with small incidence angles 
of the transmitted signal on the obstructing wall. 

III. Localization algorithms 

A. Introduction 

In this Section we develop various algorithms for two- 
dimensional localization in a UWB network composed by 
N a anchors with known positions z" = [xf, yf] £ M 2 , 
i = 0,...,N a — 1 and by a single node (dubbed mobile 
station, MS, in the following) with unknown position 9 = 
[x, y] T £ M. 2 . Any couple of the given (N a + 1) devices 
can operate in a LOS (NLOS) condition with probability P L0S 

"In the rest of the document, the subscript i has been omitted for simplicity 
when not strictly necessary. 



(1 — P LOS )- The localization algorithms described below try 
to mitigate the effects of the NLOS bias error and aim at 
generating an estimate 9 of 9 minimizing the mean square 
error (MSE) Eg < \\8 — 6\\ >. It is also important to point 
out that localization algorithms developed for NLOS scenarios 
usually consist of two steps. In fact, first the NLOS bias is 
estimated for each involved link and is used to remove the 
bias contribution in the acquired data; then the new data set 
is processed by a least-square (LS) procedure generating an 
estimate of 9 (e.g., see Q, Q, g), 0). In the following 
instead, implicit estimation of the bias for each link is adopted; 
this approach is motivated by the fact that the estimation of the 
bias for the i-th link can benefit from the information acquired 
from the other (N a — 1) links; note that in (2) and bias 
mitigation performed in a link-by-link fashion is exploited to 
assign a weight in a weighted least-square (WLS) step but 
such an approach is heuristic and differs from our ML-based 
approach. 

B. Maximum Likelihood Estimation 

If the links between the MS and the N a different anchors are 

assumed mutually independent, the ML estimation strategy of 

the unknown MS position 9, given a TOA estimate and a set of 

additional signal features for each link, can be formulated as 

9 = argmaxlnniIV 1 /r,x(ri,x i ;0), where 9 = [x, y] T £ 

e 
M 2 denotes the MS trial position, n and Xj are the TOA and 

Nf -dimensional signal vector collecting the received signal 
features acquired for the i-th link and / T ,x("; S 0) is me joint 
pdf of the TOA and the vector of features parameterized by 
the trial position 9. In the following Paragraphs the impact of 
various possible options for the ML strategy are investigated. 
1) ML estimation strategy for different sets of observed 
data: in this Paragraph two different options are considered 
for the set of features processed by the ML strategy. 

a) Option A: in this case it is assumed that the vector of 
features referring to the i-th link is x. ( = [r ma ,x,i, r m ,i, Tds,i] T > 
this choice is motivated by the large correlation between 
these random variables and the link bias b (see Paragraph 
II-B I. Then, the joint PDF f T .x{T~i, x^; 9) appearing in the ML 



strategy can be expressed as (see Eq. (fT)): 
/ T , x (ri,Xi;0) = (f b ,i «) f w ) 
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(because of the models previously described for r% lax , t™ 
and TjJ), di(9) is the distance between the i-th anchor and 
the MS trial position and fb t a (8> fw denotes the convolution 

between the joint pdf fbz(-) and the observation noise pdf 

/ u ,H = (2 7 r ( 7 2 J )- 1 /2 exp( _ w 2 /(2(T 2 i)) ^ 

Note that the shape of the function fb,z(-) under the hypoth- 
esis of LOS conditions (H LOS event) is substantially different 



from that found in NLOS conditions (H NL0S event) [4|; for this 
reason, we estimated this function in both cases applying the 
procedure described in [1] to the data collected in our mea- 
surement campaign; this led to two distinct multidimensional 
histograms, which approximate the pdf's /(,,£ (6, Xi|if L0S ) and 
/bx (6, Xi|iJ NLOS ) with a certain accuracy depending on: a) the 
quantity of acquired data; b) the sizes Ab, Ar max , Ar m and 
At^s of the quantization bins adopted in the generation of 
the histograms. Note that these sizes need to be accurately 
selected, since large bins imply a coarse approximation of 
pdf's, whereas excessively small bins require a huge amount 
of data. 

Given an estimate of the above mentioned couple of pdf's, 
the required pdf fb,ic(') can be trivially evaluated using 
the law of total probability and the P LOS — Pr{_ff LOS } and 
■Pnlos = Pr{^NLos} = 1 - Pr{i? L0S } probabilities when they 
are available, or assuming P LOS = P NL0S = 0.5 if no a priori 
information about the LOS/NLOS conditions are available. 

b) Option B: in this case the set of features employed in 
ML estimation consists of a single element, namely the delay 
spread (which exhibits the largest correlation with the NLOS 
bias; see Table lip, so that Xj = T<j s ,j and the pdf / T)X (tj, x,; 9) 
of the ML strategy becomes (see d2])) 



/t,x(7i,Xj;0) = (fb,r™ <8> f w ) n 



di{6) r ds 



co 



di{9) 



(4) 

Like in the previous case, the pdf fb. T m (') has to be estimated 
in the LOS and NLOS scenarios (see Fig. [2]i from the data 
acquired in our measurement campaign. Note that this option 
leads to a ML localization algorithm which is substantially 
simpler than that proposed in the analysis of option A. 

2) Parameterization of the observations: the ML strategies 
developed above are based on the joint pdf's /& iT .o . r ™. T ™ (•) 
and / bT ™(-) which refer to a set of distance-independent 
parameters. Since the parameters r® max , r™ and t™ s exhibit 
a lower correlation with b than their distance-dependent coun- 
terparts r max , T m , Td s , we believe that the use of the joint 
PDF 



fr,x( T ii X *; 0) = (fb,x <s> /« 



c 



(5) 



in place of d2ji deserves to be investigated (similar comments 
hold for Q). 

3) Estimation of joint pdf 's: as already explained above, the 
joint pdf's involved in the proposed ML localization strategies 
can be easily estimated from the acquired data using a simple 
procedure based on dividing the space of observed data in 
a set of bins of proper size. Such a procedure generates an 
histogram, which, unluckily, entails poor localization perfor- 
mance if employed as it is, because of the relatively small 
number of bins (adopted to avoid empty bins). To mitigate 
this problem either interpolation followed by low-pass filtering 
can be applied to raw experimental histograms (this leads 
to interpolated-histogram estimators) or the raw data can be 
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Compute the posterior probabilities: 
Pl = 5 f„Jb.x\H ws )db 
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Recompute the bias estimate: 
b = argmax/j x (A,x \H = H) 



Figure 3. Flow diagram of the proposed iterative estimator. 



fitted with analytical functions (this results in fitted-histogram 
estimators). 

C. Iterative Estimation 

Recently, an iterative estimator of both the channel state 
(i.e., LOS or NLOS conditions) and the NLOS bias b has 
been proposed in |4|. In Fig.[3]the flow diagram of a modified 
version of this iterative algorithm, employing the joint pdf's 
described above and extracted from our experimental database, 
is proposed. Note that this algorithm operates in a link-by-link 
fashion (in the diagram of Fig. [3] the link index % has been 
omitted to ease the reading). 

IV. Numerical results 

Our experimental database has also been exploited to as- 
sess the RMSE performance of the proposed algorithms for 
localization and NLOS bias mitigation schemes via computer 
simulations. In our simulations the MS coordinates are always 
(0;0); then, following (2, for the i-th anchor a received 
waveform from either the LOS database (with probability 
^i.os) or from the NLOS database (with probability P NLO s) 
was drawn randomly and was associated with the position 
zf = I di sin(27r^i); di cos(27r^i) J, where di is the dis- 
tance measured for the selected waveform. Note that these 
waveforms already include the experimental noise and thus 
no simulated noise was imposed on the waveforms (so that 
the signal-to-noise ratio is the experimental one). Finally, 
the parameter N a has been set to 3 (worst case which still 
theoretically allows unambiguous localization). 

The RMSE performance of the following algorithms has 
been evaluated: 

1) LS - A standard LS estimator for LOS environments; 
the estimation strategy can be expressed as 9 = 
argmirig 2~2i=o (c-oTj - di(9)) 2 . 

2) VE - A LS estimator exploiting TOA measurements 
corrected by the algorithm proposed by Venkatesh and 
Buehrer in Q; this algorithm relies on a statistical 



modeling of the propagation environment based on our 
experimental database. 

3) ML-4D - A ML estimator based on (J2|i and employing 
an interpolated histogram, in the likelihood function, 
referring to a distance-dependent parameterization. 

4) ML-2D - A ML estimator based on (Hjl and employing 
an interpolated histogram, in the likelihood function, 
referring to a distance-dependent parameterization. 

5) ML-2D-ID - A ML estimator based on Q and employ- 
ing an interpolated, histogram in the likelihood function, 
referring to a distance-independent parameterization. 

6) ML-4D-F - A ML estimator based on p) and employing 
a fitted histogram, in the likelihood function, referring 
to a distance-dependent parameterization. 

7) ML-4D-IT - A LS estimator based on §2§ and exploiting 
TOA measurements corrected by the modified iterative 



8) 



algorithm illustrated in Section III-C| 
ML-2D-IT 



A LS estimator based on (HI) and exploiting 
TOA measurements corrected by the modified iterative 
algorithm illustrated in Section |11I-C| 

In estimating the RMSE performance of the ML algorithms 
listed above the likelihood functions were always evaluated at 
the vertices of a square grid characterized by a step size equal 
to 10 mm. Some numerical results are compared in Fig. HI 
which illustrates the RMSE performance versus the probability 
P LOS . These results evidence that: 

1) The simple LS algorithm is outperformed by all the other 
algorithms when P LOS < 0.9; this is due to the fact that 
this strategy does not try to mitigate NLOS bias. 

2) The VE algorithm performs well at the cost of a rea- 
sonable complexity, but offers limited bias mitigation 
when P L0S = 0; in this case the ML-2D, ML-2D-ID and 
ML-4D-F estimators perform much better. 

3) The exploitation of a large set of received signal features 
does not necessarily allow to achieve better accuracy 
than a subset of them (see the curves referring to ML- 
2D and ML-4D estimators); this is due to the fact that 
the correlation between the different couples of extracted 
features is typically large, so that they provide strongly 
correlated information about the NLOS bias. 

4) Distance-dependent parameterization provides better ac- 
curacy (see the curves referring to the ML -2D and ML- 
2D-ID estimators); this can be related to the fact that r^™ 
is less correlated with the NLOS bias than its distance- 
dependent counterpart r__. 

5) The ML-4D-F estimator performs better than the ML-4D 
estimator in NLOS conditions; this means that the use of 
fitted histograms entails an improvement of localization 
accuracy. 

V. Conclusions 

In this paper various UWB localization techniques process- 
ing multiple features extracted from the received signal to mit- 
igate the problem of NLOS bias have been described and their 
accuracy has been assessed exploiting the experimental data 




Figure 4. RMSE performance versus P L os offered by various localization 
algorithms. 



acquired in an measurement campaign. Our results evidence 
that: a) a restricted set of features has to be employed; b) the 
use of distance-dependent features and of fitted histograms 
provides better performance than that offered by distance- 
independent features and interpolated histograms. 
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